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Experimental evidence of a natural parity state in ^^Mg and its impact to the 

production of neutrons for the s-process. 
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We have studied natural parity states in '^^Mg via the ■^■^Ne(®Li,d)'^^Mg reaction. Our method 
significantly improves the energy resolution of previous experiments and, as a result, we report the 
observation of a natural parity state in ^^Mg. Possible spin-parity assignments are suggested on the 
basis of published 7-ray decay experiments. The stellar rate of the ^^Ne(Q,7)^®Mg reaction is reduced 
and may give rise to an increase in the production of s-process neutrons via the ^^Ne(Q,n)^^Mg 
reaction. 

PACS numbers: 



I. INTRODUCTION 

The ^^Ne(Q:,n)^^Mg reaction is regarded as the main 
neutron source for the s-process in core He-burning mas- 
sive stars []J and is of relevance in He-shell burning in 
AGB stars [J . Our current understanding of its rate is 
one of the most important sources of uncertainty in the 
nucleosynthesis of heavy elements. 

The species ^^Ne is produced in helium-rich envi- 
ronments from ^''N, a product of the CNO cycle, via 
"N(a,7)i8F(/5+i/)i80(a,7)22Ne. The production of neu- 
trons by a-particle capture on ^^Ne occurs through a res- 
onant process involving the formation of the ^^Mg com- 
pound nucleus in an excitation range of high level density. 
The populated resonant levels decay by neutron emis- 
sion to the ground state of ^'^Mg. Figure [T] shows both 
the formation channel (open above i?a; = 10,614.78±0.03 
keV) and the outgoing neutron channel (open above 
£::r = ll,093.07±0.04keV). 

The 22Ne(Q;,7)^^Mg reaction competes with the 
^^Ne(Q;,n)2^Mg process Q above the neutron threshold. 
At temperatures of relevance to the main and weak com- 
ponents of the s-process, the neutron yield is defined by 
both the abundance of ^^Ne and the branching ratio be- 
tween the competing 7-ray and neutron exit channels 
Direct measurements of both processes have been re- 
ported by severalgroups for E'-^'^ > 800 keV (For ex- 
ample, see d, H, UH, Q and other references therein). 
However, due to the Coulomb barrier, the cross section 
remains uncertain for lower energies. A direct measure- 
ment for EI^'' < 800 keV still needs to be done. 

There is only scarce information on levels of ^^Mg 
between the a-particlc and neutron thresholds. Giesen 
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FIG. 1: Level scheme (not drawn to scale) of ^^Mg showing 
the ■^■^Ne + a entrance channel and the two competing exit 
channels ^^Mg -|- n and ^^Mg -I- 7. The levels of interest 
to this work (all below the neutron threshold) are shown. 
Note the negative Q-value (Q = -478.3±0.04 keV) for the 
^^Ne(a,n)^'^Mg process. 



et al. [vj report two natural parity resonances (at 
i;:r = 10.694(20) MeV and ^^^ = 10.949(25) MeV) in their 
22Ne(^Li,d)^^Mg experiments. As discussed by Karakas 
et al. |10|, the main source of uncertainty in the rate for 
the 22Ne(a,7)^^Mg reaction at these energies results from 
the unknown spin of the ^^Mg state at iJ^^lO. 949(25) 
MeV. However, in their calculation of the rate, Karakas 
et al. have considered only two states while more than 
20 levels are listed in the compilation of Endt [ll| for 
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excitation energies between the a-particle and neutron 
emission thresholds. A detailed study of levels in ^^Mg 
and their spins and parities at these energies is urgently 
needed. Here we report the first experimental step to- 
wards a complete understanding of the reaction rate at 
stellar temperatures. 

Both ^^Ne and ^He have a ground state with J'^=0"*'. 
Thus, preferentially natural parity states in ^^Mg can 
be populated via the "^He+^^Ne process [Hi- We stud- 
ied the astrophysically relevant natural parity states in 
26Mg between E^ = 10615 keV and E^=11093 keV via the 
^^Ne(^Li,d)^^Mg reaction, which preferentially populates 
natural parity states for direct a-particle transfer [l3| . 



II. EXPERIMENT 

Five ^^Ne targets were prepared by implanting ^^Ne 
into 40 ^g/cm^, 99.9% ^^C-enriched foils. The foils were 
floated from glass slides in deionized water and mounted 
onto aluminum frames. It is thought that the mechanical 
stability of thin foils can be improved by exposing them 
to an intense burst of light [14]. We flashed the unim- 
planted foils with a photographic strobe and observed 
their slackening. The ^^Ne beam was produced by the 
200 keV Eaton Ion Implanter at the University of North 
Carolina; two energies (20 keV and 35 keV) were used to 
implant both sides of the foils. ^^Ne was implanted in this 
way to achieve a total density of 19-22x10^^ atoms/cm^ 
(~ 7 /ig/cm^). The dose of implanted ions was estimated 
by integrating the beam current at the target, which to- 
gether with the beam stop and the implantation chamber, 
acted as a Faraday cup. Secondary electrons were sup- 
pressed with a negative voltage applied to a copper pipe 
placed in front of the target and coaxial to the beam. 
The copper pipe was cooled with liquid nitrogen to pre- 
vent natural carbon build up on the targets. It was im- 
portant to keep ^^C contamination of the targets to a 
minimum because deuterons from the ^'^C(®Li,d)^^0 re- 
action would have posed a major source of background 
in our ^^Ne(^Li,d)^^Mg spectra. An additional source of 
background could have been ^°Ne from the ion source at 
the implanter. The natural neon gas used to produce the 
beam was mass analyzed using a magnet with a mass 
resolution better than AM/M=0.01 allowed us to get an 
excellent separation of ^^Ne and ^°Ne. The ^^Ne beam 
current was kept below 400 nA to minimize heating of 
the carbon foils. 

Both the ^^Ne(^Li,d)^^Mg experiment and analysis of 
the target composition were performed with a ^Li beam 
produced by the Wright Nuclear Structure Laboratory 
ESTU Tandem van de Graaff accelerator at Yale Univer- 
sity. The reaction products were momentum-separated 
with the Enge split-pole spectrometer and detected at 
the focal plane with a position-sensitive gas ionization 
chamber and a scintillator [Tsj . 

The detector setup allowed us to separate different par- 
ticle groups by means of the magnetic rigidity (Bp), the 



energy loss (AE) measured from the cathode and the 
residual energy (E) deposited at the scintillator. The 
magnetic rigidity (Bp) was derived from independent po- 
sition signals along two parallel wires (front wire and rear 
wire). The horizontal acceptance was ±10.7 mrad, which 
corresponded to the minimum aperture available. 

The target content analysis was performed with a 7.7 
MeV ^Li beam. To calibrate the elastic scattering po- 
sition spectrum as a function of mass, we used a target 
consisting of a Si02 layer deposited on a 40/Ltg/cm^ natu- 
ral carbon substrate. Elastic scattered ^Li off the target 
was observed in three major groups, each corresponding 
to ^^C, ^""O, and ^^Si. A ^^Ne-implanted and a non- 
implanted ^^C foil were exposed to a ^Li beam as well 
(see figure 

The (®Li,d) experiment was performed with a 30 MeV 
^Li beam. We took spectra with the non-implanted tar- 
get for magnetic rigidity calibration purposes and for 
comparison with the ^^Ne- implanted target. Because the 
cross section is expected to be larger at small scatter- 
ing angles, we placed the spectrometer at 6° in the lab- 
oratory. The spectrometer field was set to center the 
deuterons from the ^^C(^Li,d)^^0 reaction at the posi- 
tion spectrum. 

We acquired deuteron position spectra with an average 
beam current of 80 pnA with the ^^Ne target and a total 
solid angle of 1.5 msr. The result is shown in figure [3l 
The energy region of interest for deuteron groups corre- 
sponding to ^^Mg states is located between the two ^^O 
doublets. Deuteron groups unobserved with the ^^C tar- 
get appear both inside and outside the region of interest. 

III. ANALYSIS AND RESULTS 

Although the detector was positioned so that the focal 
plane coincided nominally with the front wire, the two 
independent position measurements allowed us to more 
accurately locate the true focal plane by reconstructing 
the particle trajectories p^ . thereby improving the posi- 
tion resolution. Let Pi and P2 be the positions measured 
at the front and rear wire, respectively, and S the dis- 
tance between the wires {S — 2.5 cm). The trajectory 
(x, y) of particles traveling in a plane parallel to the two 
wires is described by the relation 

(Pi - P2)y + Sx - P2S = 0. (1) 

On the other hand, the {x, y) equation for the focal plane 
is 



(1 4- cot^ a)2 (l + tan^a)5 

where a is the angle between the focal plane and the front 
wire, and H is the distance from the focal plane to the 
origin (0,0). Solving {x,y) simultaneously for the two 
equations one gets the position of the particles at the fo- 
cal plane (see figure S]). The best position resolution was 
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FIG. 2: Elastic scattering off a non-implanted ^ C substrate (upper panel) and off a Ne-implanted ^ C substrate (lower 
panel). The arrow shows the position where contaminating ^"Ne would have appeared. 



obtained for S/H=2. The resulting spectrum is shown in 
figure [3l[c) ; here the energy resolution improved from 80 
keV to 58 keV , as measured for the ^^O peak at Ej;=7.12 
MeV. 

The S/H value was chosen to optimize the resolution 
of the ^^O peaks in the position spectrum. The main 
source of background was the ^^C substrate in the target, 
as can be noted by comparing panels (a) and (c) in figure 
[31 Therefore, the best peak to background ratio for the 
^^Mg states was obtained by focusing on the ^^O peaks 
instead of the ^^Mg states themselves. Nevertheless, the 
width of ^®Mg peaks for S/H=2 was still within 5% of the 
optimized value for the focused ^^Mg peaks (S/H=1.67). 
The width observed for the ^^Mg peaks is listed in table 

m 

The origin of the deuteron groups observed with the 
^^Ne target was established on the basis of a target con- 
tent analysis via ^Li elastic scattering. The target con- 
tent analysis was performed by moving the spectrometer 
to 20° in the laboratory frame and reducing the mag- 



netic field to 7.7 kG, so that elastically scattered ^Li 
ions were centered at the front position spectrum. The 
results for both the non-implanted and ^^Ne-implanted 
^^C substrates are shown in figure [21 The major peak 
appearing in the spectra corresponds to the ground state 
in ^^C. Other elastic scattering groups were also identi- 
fied. A comparison between these two elastic scattering 
spectra shows that the only group observed in the ^^Ne- 
implanted target and not observed in the non-implanted 
target is that corresponding to ^^Ne. The position where 
^"^Ne contamination from the implantation process would 
have appeared is marked as well. From figurc[21it is clear 
that ^^Al is also increased relative to the ^^C group be- 
tween both targets (probably as a result of the implan- 
tation process). Therefore, we compared the deuteron 
spectrum obtained from an ^^Al target with that of 
the ^^Ne-implanted target. We scaled the background 
with the ^^Al content determined via elastic scattering 
in the ^^Ne-implanted target (see figure [51 for a content- 
normalized comparison between spectra from both tar- 
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FIG. 3: Deuteron position spectra for 6iab~Q°- (a) The top panel shows the deuterons from the ^^C substrate. Six groups are 
observed and correspond, from left to right, to the = 8.8719(5) MeV state in ^''O (from ^2C(®Li,d)^*'0), the = 5.6214(17) 
MeV state in ^''Ne (from i''0('^Li,d)2°Ne), the £^=7.11685(14) and 6.9171(6) MeV doublet in ^''O, and the £^=6.12989(4) and 
6.0494(10) MeV doublet in ^"^O, both from ^^C(^Li,d)^''0. (b) The middle panel shows the deuteron front position spectrum 
from a ^^Ne-implanted ^^C substrate. The energy region of interest in this work is the window located between the two ^^O 
doublets. Deuteron groups not observed with the ^'^C target are seen here, (c) The lower panel shows the deuteron spectrum 
after reconstructing the focal plane for the ^^Ne- implanted ^^C substrate (see section |lIT|. The energies of some of the observed 
^®Mg states are shown. The two vertical lines correspond, from left to right, to the positions of the neutron and a-particle 
thresholds. The arrow labeled with an asterisk (*) corresponds to i?a; = 10.694 MeV. A state was not observed at this energy 
and angle (see text for discussion). The inset shows the Gaussian fits to the two peaks in the region between the neutron and 
a-particle thresholds. 



gets). We concluded the peaks observed between the two 
^^O doublets originated from the 22Ne(^Li,d)^'^Mg reac- 
tion, with a contribution of 5% from ^^Al to the area 
under the peaks. 

From the position spectrum, a magnetic rigidity and an 
excitation energy calibration can be obtained by assum- 
ing a polynomial relation for the radius of curvature p. 
This requires that at least three peaks be unambiguosly 
identified. Here we used six states that were populated 
via the i2c(6Li,d)i^O and i6o(6Li,d)20Ne reactions (see 
table m. 

The position of peaks observed was determined by fit- 



ting with a Gaussian template. The error bars of cen- 
troid positions were determined by sampling the area 
and width of the Gaussian determined by minimizing the 
value of for the fit, within error bars. A Monte Carlo 
sampling produced a set of centroid positions that deter- 
mined the size of error bars. Table |I] shows the centroids 
of the fitted peaks and the published excitation ener- 
gies for the six calibration peaks. We fitted p, the 
deuteron trajectory's radius of curvature, as a function 
of the front wire position Pi with an expression of the 
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TABLE I; States observed in this work 



Centroid 


Ex 




Nucleus Peak width 


(cliannel) 


(MeV) 


(MeV) 


(keV) 



1667.1(10) 
1792.4(15) 
1834.4(10) 
1903.2(50) 
1931.8(40) 
1996.8(20) 
2012.6(20) 
2173.5(80) 
2221.4(110) 



5.7877(26)" 
7.11685(14)" 
6.9171(6)" 
10.953(25) ' 
10.808(20) ' 
6.12989(4)" 
6.0494(10)" 
9.57(4) ' 
9.32(6)' 



25.450(11) 
26.122(1) 
26.345(2) 
26.719(61) 
26.872(50) 
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27.314(5) 
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28.46(18) 
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FIG. 4: Deuteron trajectory and focal plane diagram 
(adapted from [3). 
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FIG. 5: (Color figure) . Comparison between spectra obtained 
with both a ^^Ne- implanted target (black curve) and an ^'^Al 
target (red curve). The yield for the ^'^Al target was renormal- 
ized to match the relative accumulated charge on the 22 Ne- 
implanted target and the relative ^^Al content in the targets, 
as measured with elastic scattering of ^Li. Both spectra were 
taken with a beam energy of 30 MeV and the spectrometer 
placed at 6°. 



type 



p(Fi) = ao+ai(Pi-Pi[0]) 



(3) 



Here the a^'s are the parameters of the fit. We found 
that the best fit (x^/N=0.1, where N=3 is the number 
of degrees of freedom) was obtained with a polynomial 
of degree 2. Finally, the excitation energy of ^^Mg states 
was computed from p with a kinematic analysis. The 
error bars include contributions from the uncertainty in 



IV. DISCUSSION 



The recent calculation of the rate for the 
22Ne(a,7)26Mg reaction [13] includes contributions 
from two ^^Mg states below the neutron threshold. The 
largest source of uncertainty is from the i?:r = 10. 949(25) 
MeV state; a second state at i;:r = 10. 694(20) MeV 
(Ecm=0.078 MeV) has a negligible effect on the rate, 
as it is located too far below the Gamow window 
for He-burning temperatures (the window spans from 
Ec™=0.36 to 0.57 MeV at T=2.0xl08 K). Both 
states were observed by Giesen et al. Q using the 
^^Ne(^Li,d)^^Mg reaction. However, we find evidence 
that their £'a;=10. 949(25) MeV state corresponds to at 
least two states in ^^Mg. 

Four ^^Mg states were identified in this experiment (see 
table [T|. Two states fall in the region between the two 
i^O doublets (see figure^. The first has = 10.808(20) 
MeV and the second has = 10.953(25) MeV. 

Giesen et al. 0| and Giesen [3| studied the 
^^Ne(^Li,d)^^Mg reaction with a relatively poor resolu- 
tion of 120 keV in their deuteron spectra. In contrast, for 
the present experiment the energy resolution was 63 keV 
in the region of interest, as a result of using a solid tar- 
get. We observed two ^^Mg states outside of the region of 
interest. The first, at Ej;=9.32(6) MeV, corresponds to 
the Ea;=9. 404(20) MeV state of Giesen et al, while the 
second, observed at Ea;=9.57(4) MeV, is in agreement to 
their state at Ea;=9. 586(20) MeV. 

A state at i?j; = 10. 694(20) MeV was not observed in 
our experiment, which is consistent with the 9iab—7.5° 
spectrum of Giesen There he reports a cross sec- 

tion more than one order of magnitude smaller than that 
of the ii^a;=10. 949(25) MeV state at the same angle; no 
peak corresponding to the i?2; = 10. 694(20) MeV state in 
^^Mg can be identified in the spectrum. Nevertheless, 
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one single broad peak was observed at a higher excita- 
tion energy {E^ = 10.M9{25) MeV). With improved reso- 
lution we have resolved this peak into two ^^Mg states. 
It is thus likely that the spectroscopic factor for this dou- 
blet gets contributions from both states, one of them (the 
state with Ex — 10.808(20) MeV) is at an energy too low 
(Ecm=0.193 MeV) to contribute significantly to the reac- 
tion rate of ^^Ne(a,7)^^Mg in the region of astrophysical 
interest. 

We identify these two states as follows: Walkiewicz 
et al. (19] studied the secondary 7-rays from thermal 
neutron capture for ^^Mg and observed an 8996.5 keV 
transition to the first excited state (2+) of ^^Mg. They 
report a state with iJ^, =10.8059(4) MeV and Endt [il| 
assigned J'^=(0"'"-4+) to it. This is consistent with our 
state at £^^ = 10.808(20) MeV. The "^^mi^UAf^Mg re- 
action most likely populates states with natural parity, 
so we propose this state to have J'^=0+, 1~, 2+, 3~ or 
4+. 

Giesen et al. suggested the £^ = 10. 694(20) MeV state 
to have ^=1' , 8+ or 4+. Glatz et al. f20l| observed 
this same state to decay by 7-ray emission to the 7^=5"*" 
state at Ej,^?. 395(1) MeV. Most likely, this decay cor- 
responds to an Ml transition. This is in agreement with 
the J'^=4+ assignment suggested by Karakas et al. jTo| . 

An £^^ = 10.81(6) MeV ^^Mg state was reported by 
Crawley et al. j21j| as well. They measured forward- 
angle cross sections for 201 MeV proton inelastic scatter- 
ing and observed a forward-peaked angular distribution, 
thus suggesting J'' = l+. This state would not have been 
observed in the present experiment and most likely corre- 
S2onds to the £^=10. 824(3) MeV state observed by Moss 
[221 and listed in Endt T?!- 

Our i?a; = 10. 953(25) MeV state is at an energy consis- 
tent with the £;2;=10.945(3) MeV level hsted in Endt [ill- 
The state at £a;=10. 953(25) MeV was also observed by 
Giesen et al. and they suggested J^=3^ without be- 
ing able to discard the J'^=2~^ and 4+ assignments. Glatz 
et al. [l^l also observed this state to decay by 7-ray 
emission to the £;;j;=8.625(l) MeV, 9.169(1) MeV, and 
9.383(1) MeV states with branching ratios 29(4), 61(5), 
and 10(2), respectively. Endt [11] lists the three final 
states to have J^=5~, 6~, and 6+ respectively. Giesen 
et al.'s J'^ assignments are not consistent with these 7- 
ray decays. Assuming our experimental work populated 
natural parity states in ^^Mg, only J'^=5~, 6"*", or 7~ 
are allowed. The discrepancy between these results and 
Giesen et a/.'s probably comes from the fact that their 
DWBA analysis and J'^ assignment were performed for a 
peak consisting of two unresolved states. 

We calculated the contributions from the two resolved 
resonances to the (0,7) rate by taking the total differen- 
tial cross section observed by [7,] for their i?2;=10. 949(25) 
MeV state and splitting it into two parts. The ratio of 
the two contributions was taken to be equal to the ratio 
of the areas under the peaks observed in our deuteron 
spectrum (Fig. ^ after correcting for the background 
from ^^Al, as shown in figure [H] The spectroscopic fac- 



tor Sq for each contribution was calculated by assum- 
ing a combination of spin pairs and then by fitting a 
DWBA model computed with the code DWUCK4 Q 
to the individual experimental cross sections. For the 
state at E^=10.808(20) MeV, Sa = 1.9xl0"^ while at 
E^=10.953(25) MeV, Sa=2.8xl0"^ The upper limit of 
the reaction rate was evaluated by assigning the states at 
E:r=10.808(20) MeV and E:, = 10.953(25) MeV as J'^=0+ 
and J'^=5^, respectively. For the lower limit we took 
J'^=4+ and J'^=7~, respectively. 

Comparisons of our rates and the rates calculated by 
Karakas et al. ^ to NACRE ^ are shown in figure 
[6l Relative to the NACRE rates, the main effect is a 
reduction of the rate for temperatures below 0.3 GK. 
Thus an increase in the neutron production by the (a,n) 
reaction may be expected in these stellar environments. 
Nevertheless, it is necessary to finalize the spin-parity 
assignment of all ^^Mg states contributing to the (a, 7) 
rate and that are located below the neutron threshold. 
Further experimental studies are under way. 
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